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ABSTPACT 


An invest! Ration of thermal perturbations of the solar convective zone 
via changes in the mixing length parameter have been carried out, with a 
view toward understanding the possible solar radius and luminosity changes 
that have been cited in the literature. The results show that: (a) a 

single perturbation of a is probably not the cause of the solar radius 
change given in ref. 2, and (b) the parameter W = dinRo/dinL^ can not be 
characterized by a single value, as has been implied in recent work 
(refs. 2-5) 
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INTRODUCTION 


Recent observations of possible changes of the solar luminosity (ref. 1) 
and radius (ref. 2) have spurred theoreticians to try to model the physical 
mechanisms that might produce such changes. One possible mechanism involves 
thermal perturbations of the solar convection zone. Such perturbations can 
be modeled (to first order) by perturbing the mixing length parameter o(equal 
to the ratio of the mixing length to the pressure scale height), used in the 
standard mixing length theory of convection. In this paper wc present the 
results of such an analysis. Other work in this area can be found in 
references 2 to 5. 

METHOD OF ANALYSIS 


The stellar evolution code described by Endal and Sofia (ref. 6) was used 
for this investigation, utilizing the following assumptions and input data: 

(1) Spherically symmetric Sun in hydrostatic equilibrium 

(2) Latest Los Alamos opacities and equations of state 
(ref, 7) 

(3) Nuclear reaction rates of Fowler, et.al. (ref. 8) 

(.4) Standard mixing length theory of convection. _g 

(5) The use of 700 interior zones and a small (6x10 Mq) 

static envelope, to insure numerical accuracy. 

A one solar mass ZAMS model was evolved to the present age of the Sun 
(Tj, ■ 4.7 x 109 yrs.). A mf tch to the present luminosity and radius of the Sun 
was found by varying the initial hydrogen abundance, X, and the value of a. 

The adopted values were X 0.71242 and o ■■ 2.21772. A sequence was then 
calculated such that the time step was slowly decreased until AT - 1 yr. at 
Tq. a perturbation of a in the range 0.05% to 4X was then introduced, and the 
full non-linear, time dependent evolution followed either for 200 yrs. (witis 
AT ■ 1 yr.) or for 6x10^ yrs. (allowlnja the time step to Increase), at which 
time the normal evolutionary effects dominated the solution. In order to test 
the effect of varying the size of the small ti^ ' step, similar sequences were 
calculated using timesteps of 0.75, 2, and 5 yrs. No difference in the 
resulting models was seen. The final results are discussed in the next 
section. 
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'RESULTS 


The results of this analysis are sununarized in Table I. We can divide 
Che results into three areas, which are discussed below. In the followln,^ 
discussion we adopt the following definitions: 


= ^new " *old " 

6(x) = lA(x)*loO)/x^j^^ ■ percent change in x. 


(1) Solar Lumirosity As can be seen in Table I, the percentage change 
in is of the same sign as 6(a), and of a much larger value than the 
corresponding 6(Rg)(here r solar radius at optical depth 2/3). The 
relation between 6(Lg) and 6(a) is linear for 0.75 S 6(a) < 4.00, and can be 
written down as: 6(L,) * 0.76 • 6(a). The stnall change in slope at 

6(a) * 0.75 will be discussed later. The cheracterlstic time scale to recover 
Its initial value is a thermal time scale of ta 6x10^ yrs. 

(2) Solar i'adlus: In the solar photosphere (t s 2/3), 6(r) shows an 

initial :har.oa of the same sign as 6(a), followed by a subsequent relaxation 
tov.’ard the new equilibrium radius (larger for 6(a) negative, smaller for 6(a) 
positive) on a much more rapid time scale than 6 (Lq). For 6(o) " 2, the 
time scale to cross its original value is « 200 yrs. For t. > 2/3, (r) 

simi^ly shows ei' immediate relaxation toward the new equilibrium value. As 
Figure 1 shows, 6(Rg), and 6(Lg) to a lesser extent, changes its behavior with 
6(a) at 6(a) - <^.75, The reason for this can be traced to the finite inter- 
polation scheme used for the envelope boundry conditions. Thus the values 
for 6(Rg) and 6(Lg) are probably not reliable below 6(a) ■* 0.75. 

(3) W: A parameter used by investigators in the field to characterize 
the changes of radius and luminosity due to any perturbation is 

W = din R^/d£n Lg " 6(Ro)/6(L0). Figure 2 shows a plot of Wq vs 6(o) for the 
data In Table I. Here Wq ■= W for the first model after the pertuvbation is 
IntroQuced. As noted, the values of W are not reliable for 6(a) <0.75. We 
see that (a) since 6(R^) changes on a much faster time scale than 6(L(,), W 
will change with time (and will change sign, e.g. after « 200 yrs. for 
6(a) = 2/i) , and (b) as Figure 2 shows, the value of changes with 6(a). 
Hence, there really is no single value for W, a.-5 has been Implied by earlier 
invescigations . 
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CONCLUSION 


The main conclusions that can be drawn from this Investigation are 
as follows: 

(1) d(Lg) “ 0.76 6(a), while 6(R^) displays a marked non- 
linear behavior with 6(a). 

(2) There is no one value of W which generally characterizes 
thermal perturbations of the solar convective zone. 

(3) Even for 6(a) 5 , 6(Rq) « 2.10”^ while 6(L^) « 3 I 

Thus for 6(ILj ) * 2.5x10*'^ (ref. 2), the extremely large 
6(Lo) implied shows that a single perturbation of a is 
probably not the cause of the radius change, although 
further work for 6(a) < 0.75, and for a series of random 
a perturbations Is definitely indicated. 

(4) Since our results for 6(a) < 0.75 are probably incorrect, 
a value of W near 0.075 for small 6(o) is not precluded. 

Thus our results show chat by taking Into account Che full 
time-dependent, non-linear behavior of the problem, the 
entire range of W quoted in the literature (refs. 2 to 5) 
may be generated. 
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Table 1 


6a 

6R 

o 

6L 

• 

6R 

w - 

6L 

e 

.05 

3.818x10"^ 

.0379 

1.007x10"^ 

.1 

7.644x10“^ 

.0758 

1.008 

.3 

2.303x10"^ 

.2277 

1.015 

.75 

5.802x10"* 

.5701 

1.018 

.8 

5.974x10"* 

.6055 

.987 

.85 

6.146x10"* 

.6410 

.959 

.9 

6.319x10"* 

.6764 

.934 

.95 

6.493x10"* 

.7119 

.912 

1.0 

6.667xl0"* 

.7473 

.892 

1.25 

7.577x10”* 

.9256 

.819 

1.5 

8.460x10"* 

1.103 

.767 

2.0 

10.262x10"* 

1.459 

.703 

2,5 

12.119xlO"‘* 

1.815 

.668 

3.0 

14.033x10"* 

2.173 

.646 

3.5 

16.038x10"* 

2.533 

.633 

4.0 

18.047x10"* 

2.893 

.624 
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FIGiniE CAPTIONS 


Figure 1. The relation between £(Lg) and 6 (o), (solid line) and between 
6 (R 0 ) and 6 (a), (interrupted line). The data for 6 (a) < 0.73 
is unreliable (daahed line) . 

Figure 2. liic relation between Wq and 6 (o). The data for 6 (a) < 0.75 

is unreliable (dashed line). The dashed line for a > 3 is only 
Intended as a guide for the eye. 
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